The Cerium volume collapse: Results from the LDA+DMFT approach 
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The merger of density-functional theory in the local density approximation (LDA) and many-body 
dynamical mean field theory (DMFT) allows for an ab initio calculation of Ce including the inher- 
ent 4/ electronic correlations. We solve the DMFT equations by the quantum Monte Carlo (QMC) 
technique and calculate the Ce energy, spectrum, and double occupancy as a function of volume. At 
low temperatures, the correlation energy exhibits an anomalous region of negative curvature which 
drives the system towards a thermodynamic instability, i.e., the 7-to-a volume collapse, consistent 
with experiment. The connection of the energetic with the spectral evolution shows that the physi- 
cal origin of the energy anomaly and, thus, the volume collapse is the appearance of a quasiparticle 
resonance in the 4/-spectrum which is accompanied by a rapid growth in the double occupancy. 
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PACS numbers; 71.27.-)-a, 71. 20. Eh, 75.20.Hr 

Cerium exhibits the well known 7-a phase transition 
characterized by an unusually large volume change of 
15% 0. Similar volume collapse transitions are observed 
under pressure in Pr and Gd (for a recent review see 
[^). It is widely believed that these transitions arise 
from changes in the degree of 4/ electron correlation, as 
is reflected in both the Kondo volume collapse |^ and 
the Mott transition 0| models. The former ascribes the 
collapse to a rapid change in the valence-electron screen- 
ing of the local 4/-moment, which is accompanied by the 
appearance of an Abrikosov-Suhl-like quasiparticle peak 
at the Fermi level, lying between the remaining Hubbard- 
split 4/ spectral density. In other words, the Kondo tem- 
perature of the a-phase is much larger than that of the 
7-phase. While originally formulated in terms of the An- 
derson impurity model, similar rapid thermodynamic and 
spectral changes are seen for the lattice version, or peri- 
odic Anderson model ||] . 

The Mott transition model envisions a more abrupt 
change from itinerant, bonding character of the 4/- 
electrons in the a-phase to non-bonding, localized char- 
acter in the 7-phase, driven by changes in the 4/-4/ inter- 
site hybridization. Thus, as the ratio of the 4/ Coulomb 
interaction to the 4/-bandwidth increases with increasing 
volume, a Mott transition occurs to the 7-phase. While 
originally motivated by the Hubbard model, most recent 
support for this perspective has come from orbitally po- 
larized 1^ and self-interaction corrected modifications 
of local-density functional theory. It has been argued 
that these modified local-density methods resemble 
static mean field treatments in which the 7-phase is spin 
and orbitally polarized such that the 14 4/-bands are 
Hubbard split into one band below and 13 above the 
Fermi level for Ce. At smaller volumes, the mean-field 
polarization disappears and, thus, the a-phase resembles 
the ordinary LDA solution with all 14 bands grouped to- 
gether just above but slightly overlapping the Fermi level. 



While this approach gives the correct energy in the limit 
of small volume, and also at low temperature for large 
volume (though not the correct paramagnetic phase), a 
true many-body solution would allow for a central quasi- 
particle peak in the presence of the Hubbard splitting as 
observed in the Mott transition of the one-band Hubbard 
model 1^,^. The Kondo volume collapse calculations [||, 
on the other hand, take such electronic correlations into 
account but are based on simplified models. 

In this situation, the recently developed merger |l^] 
of LDA ||ll| and DMFT (l|J| offers an ideal means 
to study Ce realistically, including the critical intra- 
site 4/ electron correlations Only two very recent 

LDA-I-DMFT calculations have been reported to date for 
/-electron systems. Savrasov, Kotliar, and Abrahams 
(Q have employed an interpolation scheme to the DMFT 
self-energy inspired by the iterative perturbation theory 
(IPT), and present total energy calculations and the spec- 
trum for Pu. Zolfl et al. have used the non-crossing 
approximation (NCA) to report the first Ce a- and 7- 
phase spectra. Both papers find Hubbard splitting, with 
an additional quasiparticle peak at the Fermi level, for 
the respective low-volume a-phases of the two materi- 
als. For LaTiOa, a metal close to a Mott transition, 
LDA+DMFT calculations employing the IPT and NCA 
approximations to solve the DMFT equation have been 
compared to the more rigorous QMC treatment and 
quantitative differences have been reported JlTf . 

In order to analyze the Ce volume collapse, the present 
paper reports the first LDA-I-DMFT(QMC) calculations 
of the Ce total energy, investigating a wide range of vol- 
ume and temperature. We find the low-temperature total 
energy to exhibit a distinctive feature which is consistent 
with the observed Ce volume collapse. Additional calcu- 
lations of the spectral function, imaginary time Green's 
function, and double occupancy show that the energy fea- 
ture coincides with the rapid growth of both the quasi- 
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particle peak and the double occupancy. All these signa- 
tures moderate with increasing temperature. To provide 
insight into the modified local-density methods, we also 
contrast these LDA-I-DMFT(QMC) results with Hartree- 
Fock (HP) or static mean field solutions of the same 
Hamiltonians. 

The DMFT(QMC) calculations in this paper solve the 
Hamiltonians 

}^.hn,l' ni' ,a 

+ 2^1 ^ flifmaflifm'a' , (1) 
i,ma,m' a' 

where k are Brillouin zone vectors, i are lattice sites, Im 
denote the angular momentum, a is the spin quantum 
number, fiifma = clfma'^ifma^ ^^'^ t^e prime signifies 
ma ^ m'a'. The 16 x 16 matrices Hi^^p^(k) denote the 
matrix elements of the LDA Hamiltonian w.r.t. the 16 
orbitals (6s, 6p, 5d, 4/) for fee Ce, as described in Sec. 4.2 
of Ref. [|. The 4/ site energies are shifted to avoid dou- 
ble counting of the 4/-4/ Coulomb interaction which is 
explicitly incorporated in Eq.(^ via the Uf term. These 
shifted site energies and the screened if-if Coulomb in- 
teractions Uf were obtained by companion constrained- 
occupation calculations, and their values together with 
the effective 4/ electron bandwidth are shown in Fig. 5 
of Ref. 1^ as a function of volume. Since the 4/-orbitals 
are well localized, uncertainties in Uf and the 4/ site 
energies are relatively small, and, given the significant 
volume-dependence of the 4/ site energy, only translate 
into a possible volume shift. We have not included the 
spin-orbit interaction which has a rather small impact on 
LDA results for Ce, nor the intra-atomic exchange inter- 
action which is less relevant for Ce as occupations with 
more than one 4/-electron on the same site are rare. 

For the present paramagnetic calculations we take the 
4/ self-energy matrix to be diagonal 'E{iLu) Sma;m'a' and 
use two complementary approaches to perform the trans- 
formation from G{t) to G{iujn)- The approach described 
in jl^ is used on a volume subgrid to validate a new faster 
approach which fits the G(t) data with basis functions 
of the form e~'^'^'/(e~^'/^ -I- 1) and which is employed 
for the full temperature and volume grid. Unless noted 
otherwise, our results are all extrapolated to the limit of 
zero imaginary time discretization At ^ in the QMC. 
The DMFT energy per site was evaluated from 

i?DMFT-^^Tr(ifV(k)Gk(zo.„))e^""°^ +C//d. (2) 

nko" 

Here, Tr denotes the trace over the 16 x 16 matrices, T 
the temperature, N the number of k points, and 

^ ~ n/ , , ,{T^ifmaT^ifm'cr') (3) 



is a generalization of the one-band double occupation 
for multi-band models, calculated directly in the QMC 
and related to the local magnetic moment via (toz^) = 
E™.(^./m<.)-(2/13)d. 

Fig. 0a shows our calculated DMFT(QMC) energies 
-Edmft as a function of atomic volume at three temper- 
atures relative to the paramagnetic HF energies i?pMHF 
of Eq.(|^), i.e., the energy contribution due to electronic 
correlations. Similarly given are the polarized HF ener- 
gies which reproduce Sdmft at large volumes and low 
temperatures. With decreasing volume, however, the 
DMFT energies bend away from the polarized HF so- 
lutions. This striking effect becomes more pronounced, 
and begins at slightly larger volume, as temperature is 
decreased. At T = 0.054 eV, a region of negative curva- 
ture in Ey)mft — ^'PMHF is evident within the observed 
two phase region (arrows). 

Fig. [lib presents the calculated LDA+DMFT total en- 
ergy Etot{T)^Ei^DA{T) + EoMFT{T)-En,LDAiT) where 
-E'mLDA is the energy of an LDA-like solution of the 
model Hamiltonian in Eq. (|l|) ||l^ . Since both i?LDA and 
^-PMHF — ^-niLDA have positive curvature throughout the 
volume range considered, it is the negative curvature 
of the correlation energy in Fig. |^a which leads to the 
dramatic depression of the LDA-I-DMFT total energies 
in the range V = 26-28 A'^ for decreasing temperature, 
which contrasts to the smaller changes near = 34 A"^ in 
Fig. This trend is consistent with a double well struc- 
ture emerging at still lower temperatures (prohibitively 
expensive for QMC simulations), and with it the con- 
comitant volume collapse transition. The general shal- 
lowness of our T = 0.054 eV= 632 K isotherm in Fig. ^ 
is consistent with the ~550K critical end point for the 
a-"f transition, as is the width of this region compared 
to the indicated room temperature volumes of the a and 
7 phases. We estimate entropy corrections TS would 
alter the shape of this curve by less that 0.1 eV. Fur- 
thermore, pressure-volume results (not shown) obtained 
from smoothed fits to our 632-K energies yield 70-110% 
of the observed 9.6 A'^ change in atomic volume between 
7- and a-phases at and 5 GPa, respectively, and they 
track the experimental isotherm to within 2 and in some 
cases 1 A^ from 5 to 50 GPa. Altogether, this is reason- 
able agreement with experiment given our use of energies 
rather than free energies, the different temperatures, and 
the LDA and DMFT approximations. 

To clarify the physical origin of the sharp energy fea- 
ture, we employ the maximum entropy method to study 
the evolution of the 4/ spectral function A{uj) with 
atomic volume. At 1^ = 20 A'^, Fig. |^ shows that almost 
the entire spectral weight lies in a large quasiparticle peak 
with a center of gravity slightly above the chemical po- 
tential. This is similar to the LDA solution, however, 
a weak upper Hubbard band is also present even at this 
small volume. At the volumes 29 A'^ and 34 A'^ which ap- 
proximately bracket the a-7 transition, the spectrum has 
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FIG. 1. (a) Correlation energy -Edmft — -Bpmhf as a func- 
tion of atomic volume (symbols) and polarized HF energy 
Safhf — -BpMHF (dotted lines which, at large V, approach 
the DMFT curves for the respective temperatures); arrows: 
observed volume collapse from the a- to the 7-phase. The 
correlation energy sharply bends away from the polarized HF 
energy in the region of the transition, (b) The resultant nega- 
tive curvature leads to a growing depression of the total energy 
near V = 26-28 A'^ as temperature is decreased, consistent 
with an emerging double well at still lower tempertures and 
thus the a-7 transition. The curves at T = 0.544 eV were 
shifted downwards by —0.5 eV to match the energy range. 



a three peak structure which consists of a quasiparticle 
peak or Abriksov-Suhl resonance at the Fermi energy, in 
addition to the two Hubbard side bands at about ±3 eV. 
The quasiparticle peak is seen to dramatically shrink in 
going from F = 29 A'^ to F = 34 A^, which coincides 
with the range of negative curvature in the correlation 
energy. Finally, by = 46 A'^, the central peak has 
disappeared leaving only the lower and upper Hubbard 
bands in the spectrum. 

An alternative quantity that allows the study of the 
spectrum close to the Fermi energy is the Green func- 
tion at imaginary time t — (3/2, since in the limit of low 
temperature the density of states at the Fermi level is 
given by iV(0) = -{(3/tt)G[(3/2) With increasing 

volume, the low-temperature results in Fig. ^ show first 
an increase in the quasiparticle peak around the chemi- 
cal potential due to narrowing of the LDA 4/-bands, fol- 
lowed by a sharp drop at the a-7 transition, supporting 
the results of Fig. || 

To measure the itinerant or localized character we 
study the generalized double occupancy d of Eq. (|^). 
For fully localized spins, d takes its smallest value at the 
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FIG. 2. Evolution of the 4/ spectral function A{uj) with 
volume at T = 0.136 eV (u = corresponds to the chemical 
potential; curves are offset as indicated; At = 0.11 eV~^). 
Coinciding with the sharp anomaly in the correlation energy 
(Fig. |l|), the central quasiparticle resonance disappears. 



minimal fraction of sites having two 4/-electrons which 
is drain =max(0,ny — 1) for Uf < 2. In contrast, for a 
fully itinerant system in the uncorrelated U /=0 limit, the 
maximal value d^ax = (13/28)n^ is obtained, correspond- 
ing to {rimanm'cy') — (jif /lA)'^ . With decreasing volume, 
Fig. clearly shows a dramatic increase in the ratio 
{d— dmin) / {dmax — dmin) associated to the onset of delo- 
calization. However, the delocalization is not completed 
and d is still considerably lower than its maximal value, 
even at the lowest volumes of Fig. ^p. This reflects the 
correlated nature of the a-phase with a reduced Coulomb 
interaction energy Uf d and thus, compared to the uncor- 
related static mean field solution, a lower DMFT energy 
m Fig. |I|a. The rapid increase in double occupancy im- 
plies that the local magnetic moment shows a consid- 
erable change at the transition (the increase for small 
volumes is due to an increase in ny) which is, however, 
less pronounced than in the static mean field theories. 

In conclusion, our LDA-I-DMFT(QMC) calculations 
for Ce show an anomaly in the correlation energy leading 
to a shallowness in the total energy close to the critical 
endpoint for the a-^ transition, and suggest an emerging 
double well as temperature is further decreased, consis- 
tent with the observed transition. The 4/-spectra and 
a measure of the 4/-spectrum at the Fermi level show 
Hubbard splitting in the large-volume 7-phase, with an 
Abrikosov-Suhl-like quasiparticle peak first appearing at 
the Fermi level in the transition region, and then grow- 
ing at the expense of the Hubbard side-bands with subse- 
quent compression in the a-phase. These are character- 
istic attributes of the Kondo volume collapse picture for 
Ce p . On the other hand, we also find a rapid increase in 
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FIG. 3. (a) -(/3/7r) G(/3/2) vs. volume, which at low tem- 
perature gives the 4/ spectral weight at the chemical poten- 
tial; (b) Local magnetic moment and ratio of the double occu- 
pancy vs. volume. The latter indicates that the 7-a transition 
coincides with the onset of delocalization. 



the double occupancy at low temperature in going from 
the 7- to the a-phases, which could be interpreted as 
increased itinerancy of the 4/-electrons, a tenet of the 
Mott transition picture [Q. There may well be greater 
similarities between the two scenarios than has been ac- 
cepted, as argued recently in a comparison of many-body 
solutions of the respective Anderson and Hubbard model 
paradigms for these pictures j2l] . 

Our comparison of HF to the DMFT energies also of- 
fers insight into the modified local-density calculations 
for the Ce transition , which resemble static mean 

field treatments. Our results suggest that polarized so- 
lutions can give good low-temperature energies at large 
volume for the 7-phase, however, may be offset in energy 
in the low-volume paramagnetic a-regime. Even so, it is 
tempting to speculate from Fig. |l|a that their slopes are 
correct, which would reconcile standard LDA or general- 
ized gradient extensions p3| doing so well for the volume 
dependence in the itinerant phases, even though they can 
not capture the residual Hubbard splitting. 
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